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Abstract 

The derivative expansion method has been used to solve the semiclassical 
kinetic equations of quark-gluon plasma. The nonlinear spatial damping rate, 
the imaginary part of the wave vector, for the longitudinal secondary color 
waves in the long wavelength limit has been calculated numerically. 
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The Landau damping as collisionless damping is an important collective effect in quark- 
gluon plasma(QGP); it describes how the color field is affected by media during its traveling 
in QGP. In the frame of kinetic theory, it's shown that there is no Landau damping in 
linear or Abelian dominance approximation 0~! but there is nonlinear Landau damping S. 
However, all the previous works give only the temporal damping rate, i.e., the instabilities 
with complex frequency and real wave vector values. As in the electromagnetic plasma S, 
if the imaginary part of wave vector for the secondary waves resulted from the nonlinear 
interaction in plasmas may not be zero, there may be spatial damping also in QGP. The 
collisionless damping in QGP can be completely recognized only when the temporal and 
spatial damping are known. Up to now, neither the QGP kinetic theory nor the finite 
temperature QCD has given the spatial damping rate directly. In this letter, we will work 
from the QGP kinetic equations and derive the nonlinear spatial damping rate in QGP 
finally. 

The collisionless kinetic equations for QGP are @ 



1 , 9 
— p — — 

2 dp u 
1 „ d 



p»DJ(x,p) + -p»—{F^(x),f{x,p)} = 0; 



p»D,f(x,p) - - p ^—{F^(x),f(x,p)} = 0; (1) 



p^G(x,p) + = o 



where /, /, G are the distribution functions of quarks, antiquarks, gluons in QGP, re- 
spectively. F^v, T^v represent the mean field stress tensors in fundamental and adjoint 
representation, i.e., F^ u = F® u I a , = F^ u F a . I a and F a are the corresponding generators. 
The mean color field equation coupled with the kinetic equations is 



D tl F^(x)=f(x) (2) 



j u is the color current 

g f d 3 p 



f( x ) = -~J V v [(f(x,p) - f(x,p)) + 2iI a f abc G bc (x,p)] (3) 
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where f a bc is the structure constant of SU(N C ). 

To calculate the imaginary part of the wave vector for the color field conveniently, we will 
solve the QGP kinetic equations iteratively in momentum space. Noting that the energy 
density of the color field in QGP is much smaller than the hot energy density in high 
temperature and density condition, the mean field strength may be selected as the separation 
of these scales. The derivative expansion method is an effective tool to solve the nonlinear 
equation In this method, it's essential that not only the function but also the relevant 
derivative in equation will be expanded iteratively. To discuss the spatial damping, we will 
only expand the space derivative in momentum space; the wave vector and the relevant 
functions can be expanded as 

N 

3=0 
N 

Ai = Xy^\k ( ^k( 1 V--,k (JV) ); (4) 

N 

/ = E«V (j) (k (0) ,k«,...,kW) 

j=0 

where a is a dimensionless parameter introduced to denote the order of a small quantity. / 
and G can be expanded similarly. For convenience we work in the temporal gauge, A° a = 0. 
The relation between the color electrical field E and the color vector potential A can be 
expressed as: E\ = —dA l a (x)/dt. 

The color current at any order is expanded as 

J h(n) = -f / ^yP k {(f {n) ~ / (W) ) + ^rf} (5) 

From the kinetic equations we can see that the distribution functions will fluctuate when 
the mean field is applied to a plasma system as an external field, and the fluctuating functions 
reversely influence the mean field. This self-consistent relation between the mean field and 
the distribution functions of the plasma particles are well expressed by the expansions Eqs.(4) 
and Eq.(5). The leading terms describe the linear approximation. The important dispersion 
relation obtained in the linear approximation will be changed by the nonlinear interactions 
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of the eigenwaves, i.e., wave vector corrections k^, k( 2 \ , • • •, will be given to the 
eigenwave vector k^ ^ correspondingly. As in electromagnetic plasma, if the imaginary parts 
of the wave vector correction for eigenmodes are not vanishing, there will be spatial damping. 
To simplify the calculation, we suppose that there's only longitudinal vector potential: 

A\k) = k^ l A(k)/K^ (6) 

where K<® = |k< )|. 

Inserting the expansions Eqs.(4) and Eqs.(5) into the equations (1-2) in momentum 
space, then equating the coefficients of equal power in a in two sides of these equations, one 
can obtain a hierarchy of equations. 

The first order mean field equation is 

- uj 2 A {V>h (k) — j^ h (k) (7) 
The first order transport equation for the distribution function of quarks is 
p .k^fW(k,p) + 9 - £ (p-M 0) ){A (1) (^i),^(A;2,p)} 

ki+k2=k 

+ § E pdA ( i l) (k 1 ),k?J-d;f(°\k 2 , P )} = o (8) 

ki+k2=k 

where we have defined A;^ = (w,k^) and 

The first order kinetic equations for antiquarks and gluons are similar to Eq.(8) except the 
opposite signs of the terms related to {•••,•••} for antiquarks, / and A are replaced by G 
and A for gluons, respectively. 

Assuming background configuration is local neutral, the zero order colorless distribu- 
tion functions can be chosen as the Fermi-Dirac and Bose-Einstein equilibrium distribution 
functions, respectively: 

f(°\k,p) = f (0) (k,p) = (e pp - u + l)" 1 ; G (0) (M = (e Pp - u - l)' 1 (10) 
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where (3 = 1/T is the inverse of the temperature and U is the local flow velocity(normalized 
to UfJJ^ = 1). Using Eq.(8) and Eq.(6), we can express the first order distribution functions 
in terms of the corresponding zero order distribution functions and the first order field 
potential 



f^(k,p) = -gco '' 



fW(h,p)=gu: 



p ■ fc(°) + ip°0+ dpi K(°) ' 
p o k m dfW(p)AW(k)_ 
p ■ fc(°) + ip°0+ d^ #(°) ' 
p fc(°> dG , (°)(p)^( 1 )(A;) 



(11) 



p ■ fc(°) + ip°0+ dp 4 if (°) 

By inserting Eq.(ll) into the mean field equation Eq.(7) we can obtain the dispersion relation 
satisfied by the frequency and wave vector of the eigenwave in the first order approximation. 



e[lo 



k<°)) 



^ if (0)2 



1 - 



2AT(°) 



In 



if(°) - 



CO 



oo 







(12) 



and the solution is 



A ( 1 ) cr — o " T7" 7 
k(°) ~~ _ Ul 



7T 



CO 



k(°) 



(13) 



where -E£ (0 ) and 0£ {o) are the initial amplitude and phase of the oscillation, respectively 



'k(°) 



(2N C + N f)g 2 T 2 /18 is the plasma frequency, with Nf being the number of flavors for 



quarks. The dispersion relation Eq.(12) agrees with the leading order result using the hard 
thermal loops in finite temperature QCD and the classical nature of the hard thermal loops 
has been investigated extensively by Blaizot and Jancu, etc. I§0. As shown by U. Heinz, 



the eigenwaves satisfying Eq.(12) are always timelike, i.e. 



> l(we set U — c — 1 for 



convenience). It means that the phase velocity of eigenwaves is bigger than the velocity 
of light. So these waves can't exchange energy with plasma particles and do not undergo 
damping in the linear approximation &!. In the long- wavelength region, the dispersion 
relation Eq.(12) is reduced to the form 



co 2 = co 2 p + ^ (0) 
F 5 



2 



(14) 
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The following identities will be used to simplify the color current contributed by gluons 
in the following discussion: 

fabcfabd = N c 5 c d, ifabcGbc = tr{F a G)\ [F a , F&j = ifabcFc, 

I d tr(Fd[F a , F b }) = N c [I a , h]; I d tr(F d [F a , [F b , F c ]]) = N c [I a , [l\ 1% 

tr({F a , {F b , F c }}F d }) = 45 ad 5 cb + 25 ab 5 cd + 25 ac 5 bd + N c d ade d bce (15) 

The second order equations can be discussed in the same manner. Substituting the 
results in the expression for the second color current by using the distribution functions 
obtained from the second order kinetic equations, one can find the second order mean field 
equation is reduced to 



-u 2 e{uoM 0) )A {2 \k) 



d 3 p 



1 P ■ k(°) P ■ k^ P ■ W°> 

9 V ■ + ip<>0+ K(°) k J^ 2=k p-k 2 + ip°Q+ K [ 0) K, 



(0) 



Eq.(16) describes the three-wave processes owing to the nonlinear coupling term 
] (ki) A^ '(/c 2 ) of two eigenwaves with wave vectors kx and k 2 into a secondary wave with 



wave vector k. However, as pointed out in the Refs. |pT],^2[, the three- wave processes are 
forbidden, there will be nonlinear effects of the three-wave processes on wave vectors only 
in higher order perturbation. 

Now let's discuss the third order equations. The third order field equation is 



u 2 



^ k 1 +k 2 =k J\\ &2 



+ 9 2 E E Iwtowo) ^^)^^^),^ 1 ^^)]] =j m (k) (17) 



ki+k2=k fc 3 +fc 4 =fc 2 J^-i J^2 • fv 4 

and the third order kinetic equation for quark is 



p. k WfW(k,p)-k ( . 2) p l fV(k,p)-gp l E [A?\k x )jW{k 2 ,p)]- Pi kVfW{k,p) 

ki+k 2 =k 
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-gin E [4 2) (^),/ (1) (^)] + § Yl p-k?{A?\h),d;f^(k 2 , P )} 

ki+k2=k k±+k2=k 

+§ E p-M 0) {4 2) (^i),^/ (1) (^p)}-| E ntfh^MrtfMfap)} 

k\+k2=k k\+k2=k 

+§ E MA^^O.^/^Cfe.pJl + f E ^{4 2) (^i),M?^/ (1) (A; 2 ,p)} 

fcl+A;2=fc k\+k2=k 

+§ E M^ 2) ,Ml ) ^/ (0) (fe,p)}+f e (p-knA^th),^^^)} 

k±+k2=k ki+k2=k 

-f E VikfhAf\k 1 ),d> p f { V){k 2 ,p)}- 9 - PikfhAf\k x )^\k^p)} 

k\+k2=k k\+k2=k 

+§ E ^R (1) (^i),A:S 1 >;/ (1) (A: 2 ,p)} + | E PiR^i), k^d;f°\k 2 ,p)} 

k\+k2=k k\+k2=k 

+§ E M4 (3) (^i),^;/ (0 Hfc 2 ,p)} 

fcl+fc2=fc 

-y E E P l {[A ( l\h\Af\k A )ldl^\k 2 ,p)} 

k\+k2=k k^+k^=k\ 

-y E E M^C^^fC**)],^^^)} 

-y E E ^{[4 2) (^),4 ) (A; 4 )],^/( )(A; 2 , P )} = (18) 

fcl+fc2=fc k^+k^=k\ 

The third order kinetic equations for anti-quarks and gluons are similar except for the 
changes pointed out in Eq.(8). 

Analogously to the electromagnetic plasma theory 111 , the wave vector k can be expressed 
as k = (3 — id, where (3 = \(3\ is the phase constant and a = \a\ is damping constant, i.e., the 
amplitude of the eigenwave will exponentially decrease in the direction of a. The nonlinear 
spatial damping rate at this order is defined as 



a 



-Im\^\ = -ImK {2 \ (19) 



As the oscillations are developed from random thermal motions, we have (A^ a (k)) = 0, 
where ( ) means the average with respect to the random phase of the oscillations. So we 
obtain 

{A ^(k)A^ b *(k')) = (2n)W(k - k')S b a (A^(k)) k(Q)u) ; 
(A {1)2 (k)) kWuJ = - 2 [5(uj - c k( o))]4 ( o); 4,o) = (20) 



where V is the volume of the plasma; J k (o) characterizes the total intensity of the fluctuating 
oscillation with frequency and — uj^. In an equilibrium plasma and the long wavelength 
limit (KM = 0) there are J k (o) = AttT. The average of the product of the field potential can 
be expanded as &0 

( aw (kM^Kh)^ (k 3 ) a® (fc 4 )> = (h) (k 2 )) < a® (k 3 ) (fc 4 )> 

+ (A^(k l )A^(h))(A^(k 2 )A^(h)) 

+ (A^\k l )A {l \k 4 ))(A^\k 2 )A^\h)) (21) 



By inserting the third order distribution functions Eqs.(18) into Eq.(17), multiplying 
both sides of the equation by A^ 1 ' b (k') and performing the average of the results with respect 
to the random phase, we have 



Im 



d 3 p 



. k (or 



p 



. k (2) 



E p (2vr)3 V K(°) J (p . jfc(o) + ip o Q +y 



N ^iM i N c dGi ° ){p) 



dE f 



dE n 



Ax + A 2 ; (22) 



A x = Im\ -— 



d 3 p 



d% 



■ (0), 



X 



X 



(2tc) 3 E p J (2tt) 4 ^ 'K - \ K { 

i / p°k(°)-kS 0) 



c^i y / ^ y 

(°)J \K(°)J 



p ■ (A;(°) - k^) + ip°0+ \P ■ + ipOQ- 



k (o) . k W 



P 



k (0) 



^p-kS 0) -p o k (°)- k i ) ' 

(p • fc(°) + *p°0+) 2 

-^ lP - k ( o )/p o + k (°) - k S 0) 



+ (p-k 1 (°))- 



P ■ k S 0) + i P °o+y 



(p-k) 



-cup-kS 0) /p + k(°) -kf 
(p- fc(°) + ip°0+) 2 

1 



(o) 



'7 d/W(p) 22 rfG(°>(p) \ 



+ (p. k f») (p ' k,0,) 



p° Vp • M 0) p ■ k{0) + 



'7 « 22^GW^\ 
x 6 J d£ p 2 3 d£p 2 J 



= Im< 



d 3 p 



d^kx 

(2n)*E p J (2tt) 4 
1 



M<D 9 \ ( TV ^ + TV ^-M] 



X 



X 



X 



p ■ fc(°) + ip°0^ 



v 

6g 2 /p-k(°)\ 2 /p-k(°) x 2 



Mr, 



(o) 



p. (jfe(°) - kf ] ) +ip°0+ 



s p • kx + ip°0 + p ■ + ip°0~ 



+ 



i2^ 4 / P • k(°) \ / p • kS 0) \ / p • (k(°> - k { ; 



,(0)^ 



K{ 



(0) 



p- kx+ ip°0 + p ■ + ip°0- 



d 3 P ' ( AT df^(p>) dG^(p') \ 



| k (o) _ k C0) | 
dE' p J 



X 



X 



| k (o)_ k (o)| ) { K® ) [ K { 0) ) (^-^i) 2 e (c- Wl ,k(o)-kS 0) ; 

1 / UJ\ UJ 



p> ■ (Jfc(°) - fcj 0) ) + zp°0+ \p' ■ kf ] + zp°0+ P' ■ + ip'°0+ 
To get the numerical result of the nonlinear spatial damping rate, we perform the integrals 
in cylindrical coordinates and in the local rest frame of the plasma particle. 

At first the left hand side of Eq.(22) can be simplified by averaging over all directions of 
k^ 2 ) and selecting the direction of k^ ) as the direction of the polar axis, i.e., we have 

J £ p (2ir) 3 (Al»)) 2 (p- *•<») + tp°0+) 2 \ 'dE p J W c dE T ) 

-'^f^M< f(0 ' M+N ^) (23) 

where <fi is the angle between p and k 1 - -*. ( ) means the average over all directions of 
k (2) @§ 

Now we can see that Im(K^) can be extracted from the left hand side of Eq.(22), 
and the numerical value of Im(K^) is determined by A\ and A 2 . Before performing the 
integrals in Ai and A 2 , we analyze the mechanism of the nonlinear spatial damping. The 
important relation p°/(p ■ k' + ip°0 + ) = P [1/ (a/ — v • k'] — i7i5(uj' — v • k') is very useful 
in our discussion, where v = p/p° is the velocity of particle and P stand for the principal 
value of the function. We can see that Im(K^) is not vanishing only when at least one of 
the imaginary parts of l/(p- k^ +ip°0 + ), l/(p- kf^ + ip°0 + ) and l/(p- (k^ — kf 1 ) + ip°0 + ) 
is not zero. As is shown in the linear approximation, the color eigenwaves are timelike, i.e., 
the phase velocity of waves cann't approach the velocity of particle and the imaginary parts 
of l/(p • jfeW + ip°0+) and l/(p • fc$ 0) + ip°0 + ) are zero. The term l/(p • (JfeW - k^) + ip°0 + ) 
appearing in Ai and A 2 are very crucial, it describes that the eigenwaves with wave vectors 
k(°) and kf ■* may produce the secondary waves with wave vectors k(°) — k^ through the 
nonlinear interactions. Even if the eigenwaves are timelike, these secondary waves may be 
spacelike, i.e., their phase velocity j^p^o) ma y be smaller than the velocity of light and 
can approach the velocity of particle. So the imaginary part of 5 ((oj — uoi) — v ■ (k(°) — 
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may not be zero and the secondary waves may exchange energy with plasma particles and be 
damped by the particles. It should be emphasized that though the mechanism is analogous 
to the mechanism of nonlinear Landau damping in eletromagnetic plasmas, the nonlinear 
coupling of the waves takes place through the nonlinear relations A^(ki), A^ l \k2) etc. in 
Eq.(17), i.e., it describes the nonlinear and non-Abelian characteristics of QGP. 

The similar integral to A\ and A<i has been discussed in Ref. @], by taking only the 
leading order in g and taking the long- wavelength limit, we find 

A 1 + A 2 ~ 0.42T. (24) 

After finishing the integral of the right hand side of Eq.(23), and from Eq.(24), we obtain 
the nonlinear spatial damping rate for a pure gluon gas in the long wavelength limit 

a = -Im(K {2) ) ~ 2.52g 2 T (25) 

Now we summarize briefly. In this paper, we use the derivative expansion method to 
study the nonlinear and the non-Abelian effects given by the kinetic equations of QGP. By 
solving the equations to the third order, the numerical result of the spatial damping rate 
a for the secondary waves in the long wavelength limit is obtained. This result indicates 
that the amplitude of the secondary waves in QGP decrease exponentially, and its damping 
rate is 2.52g 2 T Neper approximately. It makes us have a concrete understanding about the 
spatial traveling of color waves in QGP. 
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